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Abstract 
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nos in the final state. Initial state radiation representing the most relevant quantum 
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PROGRAM SUMMARY 



Title of program: GENTLE/4fan 
Version: 2.0 (December 1996) 
Catalogue number: 
Program obtainable from: 

[http: / / www.ith.de / theory] or on request by e-mail from the authors 
Licensing provisions: non 
Computers: all 
Operating system: all 
Program language: FORTRAN-77 

Memory required to execute with typical data: 0.5 Mb 

No. of bits per word: 32 

No. of lines in distributed program: 8300 

Other programs called: none 

External files needed: none 

Keywords: e~^e~ annihilation, QED, higher order corrections, 
four-fermion production, WW, ZZ, Z7, 77, ZH production 
Nature of physical problem: 

Description of two-boson production in e^e~ annihilation at high energies 
with selected possible background contributions; initial state QED corrections 
(ISR); selected final state corrections; anomalous triple gauge boson couplings 
Method of solution: 

Semi- analytic; numerical integration of analytic formulae: 

(i) Born approximation: Numerical integrations over two virtualities Si, S2 plus 
optionally over the boson production angle 6; 

(ii) QED corrections in the flux function approach (FF): one additional nu- 
merical integration over virtuality s'; 

(iii) QED corrections in the structure function approach (SF): alternatively 
to (ii), two additional numerical integrations over Xi,X2 

Restrictions on complexity of the problem: 

Only selective experimental cuts are possible; the program calculates total 
cross-sections and the distribution in the boson production angle; background 
contributions are available only for selected final state configurations; no in- 
clusion of pure weak corrections 
Typical running time: 

The running time strongly depends upon the options used; one finds e.g.: 

Born cross-section on-shell: below one sec, off-shell: few sees.; 

LLA QED corrections (FF) on-shell: few sees.; off-shell: few mins.; 

LLA QED corrections (SF) on-shell: few sees.; off-shell: several mins.; 

complete ISR (FF): on-shell: few mins., off-shell: several hours; 

with background: from few sees, (no QED) to several hours (with QED); 

with anomalous couplings: from few sees, (no QED) to several mins. (with 

QED) 
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LONG WRITE-UP 



1 Introduction 

With the advent of LEP2, now running above the W pair threshold, e+e~ 
coUisions have entered their multiple heavy boson production era which will 
continue with a higher energy linear e'^e~ collider (LC). Among the most in- 
teresting processes to be studied in high energy e"'"e~ annihilation is boson pair 
production. As heavy bosons decay, their physics is only accessible through 
their decay products. Thus, at LEP2 or LC, boson pair production must be 
accessed through the experimental analysis of four-fermion final states, 

e+e- ^ {W+W-, ZZ, Z7, 77, ZH, . . .) ^ /1/2/3/4 . (LI) 



2 Physics Background 

The Standard Model (SM) four-fermion final states produced in e"'"e~ annihila- 
tion have been classified by several authors [1-4]. Using the semi-analytic 
approach, we have so far studied the production of four-fermion final states 
that contain neither identical particles nor electrons or electron neutrinos. 
These four-fermion final states are printed in boldface in tables 1 and 2. 

Table 1 

Number of Feynman diagrams for 'CC ' type final states. 
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2. 1 Off-shell boson pair production 

For semi-analytic tree level calculations (flag IBORNF=0), the four-particle 
phase space is parametrized by the boson scattering solid angle, the boson 
decay solid products' angles in the corresponding boson rest frames, and the 
two final state fermion pair invariant masses si and S2 [5,6]. After analytic 
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Table 2 

Number of Feynman diagrams for NC type final states. 
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integration of all angular variables, the total cross-section for double-resonant 
four-fermion production is given by 

(7'^es(s) = y dSiPB(Sl) y dS2PB(S2) (7o(s;Si,S2) (2.1) 

with the Breit-Wigner densities 

TT |Sj - ISi YbIMbY 

(2.2) 

of the resonant bosons. Explicit formulae for the twofold differential cross- 
sections in the double-differential approximation cro(s; Si, S2) are found in ref- 
erence [7] for W pair production, in [3] for Z pair production, and in [8] 
for associated Higgs production. We emphasize that these twofold differential 
cross-sections are given by very simple formulae. The Breit-Wigner functions 
are smoothed by mapping (see appendix A). Flag IMAP allows for the CC case to 
use the Breit-Wigner functions directly. In the latter case the program works 
slower by about a factor of two. We note that in this case the invariant mass 
distributions are accessible. With flag INCPRC one may select ZZ, 77, or their 
common production. For the CC case, the on-shell limit may be calculated with 
flag IONSHL=0. 

The realization of the various generic functions mentioned here and in the 
following sections is described in section 3.3. 

2.2 Background contributions 

Although double-resonant amplitudes are dominant in most relevant cases, 
single- and non-resonant 'background' amplitudes contribute non-negligible 



5 



to the cross-sections. Including background Feynman diagrams, one can use 
the numbers from tables 1 and 2 to introduce a transparent notation [3]: 

CC Processes with final states of type /"/f/2"/2 are called charged current 
processes: CCll, CCIO, CC09; CC20, CC18. 

NC Processes with final states of type /1/1/2/2 are called neutral current pro- 
cesses: NC32, NC24, NCIO, NC06; NC48, NC20, NC21, NC19; NC4-16, NC4-12, 
NC4-03, NC12; NC4-36, NC4-09. 

mix Processes which may be considered as both CC and NC types are called 
mixed processes: inix43, mixl9; mix56. An example for a mixed process 
is the production of uddu = uudd. 

In the NC and mix classes, also Feynman diagrams with Higgs bosons con- 
tribute. In addition to the two-boson production classes CC03 [7] {W pair pro- 
duction), NC02 [6] {Z pair production), and NCOS [6] {ZZ, ^7,77 production), 
GENTLE/4fan is able to perform calculations for CCll, CCIO, and CC09 [5] as 
well as NC32, NC24, NCIO, NC06 [9]. This is chosen by fiag IPROC. In the above 
NC cases, also Higgs diagrams are included [8]. 

Complete tree level four-fermion production cross-sections are given by the 
generic formula 

Born/ N /■ 1 /■ 1 dVk(s;Si,S2) , . 

with A = A(s; Si, S2), A(a, b, c) = a'^ + b'^ + c^ — 2ab — 2ac—2bc and 

- — = Ck{s; si, S2)gk{s; si, S2) . (2.4) 
dsi dS2 

The coefficient functions Ck represent coupling constants and boson propaga- 
tors, while Qk is a kinematical function obtained after analytical performance 
of angular integration. The index k labels cross-section contributions with 
both different coupling structure and different Feynman topology. Noticing 
that often different twofold differential cross-sections d^ (7k/ dsi ds2 are related 
by symmetries in their arguments, one realizes that equation (2.3) enables ef- 
ficient coding of the cross- sect ion. For CC processes, background contributions 
are switched off with flag IBCKGR=0. With flag ICHNNL one may choose the 
final state topology in the CC case and with IFERMl and IFERM2 in the NC 
case. 
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2.3 Initial state QED radiative corrections 



In e'^e~ annihilation, the most relevant radiative corrections are initial state 
QED corrections (ISR). The phase space is now parameterized by the center- 

of-mass photon scattering solid angle, by the boson scattering solid angle in 
the two-boson (or, equivalently, four-fermion) rest frame, by the boson decay 
products' solid angles in the corresponding boson rest frames, by the final 
state fermion pair invariant masses, and by the four-fermion (or reduced) 
invariant mass s' [6]. After analytic integration over all angular degrees of 
freedom, the 0{a) ISR-corrected total four-fermion production cross-section 
with soft photon exponentiation is determined in the flux function approach 
(flag ICONVL=0) by 



with 



dVk(s,S^;Si,S2) _ , -ry 

as I dS2 as' l 



(2.6) 



where (3^ = 2 (a/vr) (L — 1), L = hi{s/ml). The soft -|- virtual and hard cor- 
rections Sk and Tik separate into a universal, factorizing, process-independent 
and a non-universal, non-factorizing, process-dependent part. The latter is 
included with flag ITN0NU=1. They are given by [6,10] 

5fe = [1 + S{s)\ ak,o{s'; si, S2) + (J^ fe(s'; si, S2) , 

Hk ^ H{s,s') ak,o{s';si,S2) cr^ ^^.(s, s'; si, S2) (2.7) 
^ ^ ' ^ ' 

Universal Part Non-universal Part 

with ak,o{s'; Si, S2) = [\/A/ ("^ts'^) GkW', Si, S2). The universal 0{a) soft-|-virtual 
and hard radiators S and H read 



(2.8) 
(2.9) 

Flag IZERO influences S{s). If the index k labels pure s-channel e'^e" annihila- 
tion, the non-universal contribution vanishes. Non-universal ISR contributions 
are analytically very complex, but, in contrast to the universal corrections, not 
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logarithmically enhanced. 



For the CC03 process, a significant QED radiative correction is the so-called 
Coulomb singularity [11,12,14] originating from the exchange of photons be- 
tween the two resonant W bosons. It effects in a factor (1 -|- J^couiomb) to be 
multiplied to the twofold differential cross-section of equation (2.4) or to the 
threefold differential cross-section of equation (2.6). Details are defined with 
fiag ICOLMB. 

The leading logarithmic ISR contributions (essentially what we called the uni- 
versal ISR contributions) may also be determined using the so-called structure 
function approach [15] (fiag IC0NVL=1). The basic formula, adapted for our 



dV, 



e+c~— >4f n-y 

dsi ds2 



1 1 

= J J dx_ D(x+)D(x_)(7o(s';si,S2) (2.10) 



with the tree level twofold differential cross-section ctq evaluated at s' = sx+X-. 
The variables x+ and a;_ represent the momentum fractions of the e"*" and e~ 
after the radiation of photons. Their minimum values are given by 



x^ — 

s 



x-^ = (y/^^+v^y. (2.11) 

x+s 

An exhaustive discussion of the application of the structure function approach 
to 4/ production may be found in [1,15] and in references quoted therein. 
Our implementation of the structure functions D{x) follows equations (4.18) 
to (4.21) of reference [5]. Note that, to improve the calculation, cto can be 
replaced by a "dressed" tree level cross-section which includes the running 
QED coupling and/or the Coulomb factor. Since the soft-photon pole is pro- 
portional to {L — 1) rather than to the leading logarithm L alone there are the 
so-called BETA and ETA alternatives for the implementation of the structure 
functions [15]. They are selected by fiag IZETTA. 



For both descriptions of QED corrections, the details of higher order terms 
are defined by fiag IQEDHS. 

In addition to the total cross-section, GENTLE/4fan can compute the angu- 
lar distribution, bin-integrated cross-sections, and moments for several other 
physically relevant quantities. These quantities are selected with the variables 
IOCS and IREGIM, the latter being chosen with fiags IRMAX, IRSTP: 
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(1) the radiative final state four-fermion invariant mass loss X = {s—s') / {2^/s) 
for the ceil family in the flux function as well as the structure function 
approach; 

(2) the true radiative energy loss X — \fsj1 [(1 — a;+) + (1 — for the 
ceil family in the structure function approach; 

(3) the W mass shift X = [{^y/s^ + ^) /2 - Mw] for the CCll family in 
the flux function as well as the structure function approach: 

(4) moments X related to the angular distribution are treated in section 2.4. 

The running of loop variable IMOMN=n is defined by flags IMMIN, IMMAX. It 

is used for the selection of the moment of n*'* degree M(x,n) of a quantity X 
(calculated with QED corrections): 

MiX,n) = ^j-^ (2.12) 

where X is selected by IREGIM and the integration is over the phase space 
variables upon which X depends and ax is the cross-section differential in 
exactly those phase space variables. Note that M(x,i) is an average of the 
quantity X. 

Finally, wc mention that an inclusive implementation of QCD corrections is 
provided. Firstly, QCD corrections to the total W width have been imple- 
mented by a factor (1 -|- 20:5/377) to be multiplied to F^^^. Secondly, the final 
state radiative QCD corrections to the formation of a quark- ant i- quark pair 
are built in by multiplying the relevant branching ratio with (1 -|- cks/tt). 



2.4 Angular distributions 



Because of the importance of the CC03 process, GENTLE/4fan is also able 
to compute W production angular distributions, bin-integrated cross-sections 
and moments for this process. A selection is possible by the flags IDCS and 
IREGIM. In the first and last of the three mentioned cases, the last analytical 
integration is not carried out. Instead, in Born approximation the analytic 
threefold differential cross-section 

, ^'^7 , = ^,j:cror{s,s2,s,;cos0) (2.13) 

dsi ds2 d cos 9 ns^ 

is used. The angular kinematic functions ^^'^°^(s, S2i Si; cos ^) were published 
in [5]. The cross-section dfJccos/ dcos^^ is obtained upon integration over Si 
and S2. In GENTLE/4f an, as a measure for the CC03 angular distribution, the 
first four moments of the differential cross-section (2.13) may be calculated as 
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its convolution with Chebyshev polynomials T„(cos^). The Born moment of 
^th is given by 



Mfcose,n) = / dsi ds2 dcosg Tn(cos6') ^ 
7 dsi ds2 d cosb* 

The zeroth moment is just the total cross-section. 



(2.14) 



Initial state QED corrections to the angular distributions and moments have 
to be calculated in the structure function approach. With photon emission, 
the W production angle in the laboratory system, 6'iab, and in the center- 
of-mass system of the four final state fermions, 9, are different. The Lorentz 
boost relating them depends on both x+ and x_ . The QED corrected angular 
moments with ISR are given by 



M, 



(cosfliab,") 



dsi J ds2 J dx+ J dx_ J d cos^ 



xD{x+)D{x^) r„[cos^iab(cos^)] 



d'^o'ccos 



dsi ds2 d cos^ 



(2.15) 



With flags IDCS, IBIN, IBINNU, the user may select the calculation of the 
angular distribution: 



d(7cc03 
d COs6'lab 



J dsi J ds2 J dx+ J dx_ 



d cos9i 



d^CTccos 



dsi ds2 d cos^i 



(2.16) 



with the threefold differential tree level cross-section d^fJccos/ ( dsi ds2 dcos^i) 
depending on s' — sx+X-, Si, S2, and cos^j(s, Si, S2j cos^iab)- Index i 

indicates that there may be zero, one, or two solutions for the Lorentz boost. 



In analogy to (2.15), the bin-integrated cross-section formulae may be written 
with the boosted angular argument under the integral, thus allowing for rel- 
atively simple analytic angular integrations (and avoiding the appearance of 
the Jacobean) if, of course, the integration limits set in the laboratory frame 
are transformed into angular limits in the boosted frame. In this respect, the 
approach to the bin-integrated cross-section differs from that for the angular 
distribution where we integrate over the laboratory angle (and have to take 
into account the Jacobean). A choice can be made with flag IBIN. 
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2.5 Anomalous couplings 

Angular distributions as defined in the above section 2.4 are especially useful 
in the context of anomalous triple boson couphngs [16]. In GENTLE/4f cin, only 

those anomalous 'yW~^W~ and ZW^W~ couplings are implemented which 
obey Lorentz invariance, conserve CP, and do not modify the electromagnetic 
interaction. These constraints allow a restricted set of anomalous contributions 
to the triple gauge boson vertex. These couplings are denoted by g^{V) with 
V e {7, Z}. The a — s,x,y,z denotes the coupling type including the SM 
triple-boson coupling g^iV). The couplings are derived from the anomalous 
couplings defined in equations (B.3) to (B.8). The threefold differential CC03 
cross-section is then given by [17] 



where and k run over s,x,y,z. The summation over s-channel photon 
and Z amplitudes is contained in the coefficient functions C. Compared to 
the SM where only three kinematical functions describe the CC03 process, the 
inclusion of anomalous couplings adds considerable complexity to the calcula- 
tion. The user may select a treatment of anomalous couplings with flag lANO. 
The potential influence of an extra heavy neutral boson Z' is explained in 
appendix C. 



3 Description of the Program 

GENTLE/4f an is a Fortran77 program for the computation of total and dif- 
ferential cross-sections for four-fcrmion production in e^e~ annihilation with 
inclusion of the relevant initial state QED corrections. It originated from older 
versions [18]. 



3.1 Features of GENTLE/4fan 

The cross-sections, angular distributions, and moments presented in equa- 



tions (2.1), (2.3), (2.5), (2.10), (2.12), (2.14), (2.15), (2.16), and (2.17) are 



computed in GENTLE/4f an by numerically integrating analytical formulae for 
expressions such as equation (2.6). The analytical input formulae are found 
in [3,5,7-9,6,10-12,17] and in references therein. Numerical integrations are 
performed by a fast and stable self-adaptive Simpson algorithm realized in the 




(2.17) 
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subroutines SIMPS and FDSIMP0. 

For transparency - and also for historical reasons - version 2.0 of GENTLE/4f an 
is built from three pieces of code. The corresponding ramification is controlled 
by the flag IPROC. The first part contains the charged current (CC) reactions 
including 0{a) initial state QED corrections and anomalous couplings. The 
second part is the 4f an code with NC reactions. The third part, steered by the 
subroutine NCGENT contains the code for the NC02 and NCOS processes with 
complete 0{a) initial state QED corrections. Subsequently we will call the 
first part the CC part, the second one will be named the 4f an branch, and the 
third part will be called NCqed part. The first and second part together are 
referred to as the CC/NC branch. 

At present, the GENTLE/4f an may treat the following four-fermion final states: 

(1) CC03 (with complete ISR, angular distributions, and anomalous cou- 
plings) [6,10,17] 

(2) NC02, NCOS (with complete ISR) [6] 

(3) CC09, CCIO, ecu [5] 

(4) NC06, NCIO, NC24, NC32 [9] 

(5) NC of (4) + Higgs [8] 

Initial state radiation is implemented. For total cross-sections universal ini- 
tial state radiative corrections as described in equations (2.5) to (2.7) are 
computed. In addition, non-universal ISR is available for the CC03, NC02, and 
NCOS processes. For the CCll family, leading logarithmic initial state QED cor- 
rections can also be calculated in the structure function approach as indicated 
in eq. (2.15). Further, the Coulomb correction to CC03 and inclusive QCD 
corrections (to both the denominators of the Breit-Wigner functions and final 
states with quarks) have been implemented. The user can choose between the 
different prescriptions of references [11,12] for the Coulomb correction via the 
flag ICOLMB. 

Cuts may be applied on the invariant masses si and S2ln the structure function 
approach, cuts on the electron and positron momentum fractions after initial 
state radiation may be imposed. In the NCqed branch of the code, a minimum 
value for the fraction s'/s may be set as well as for si/s and S2/S. A cut on 
the CC03 W production angle is not directly supported by the code but may 
be imposed using the option for bin-integration. 

In general, final state masses are neglected in GENTLE/4fan, i.e. the program 
uses the small mass approximation. Where needed, however, masses are re- 
tained in the phase space factors. Masses of heavy particles coupled to the 

^ The Fortran code FDSIMP by Yu. Sedykh is an extension of the SIMPS code by 
I.N Silin. 
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Higgs boson are taken into account where appropriate. In the NCqed branch 
of the program, masses may be chosen by altering the values of AMI or AM2 in 
the initialization routine NCINOO. 

3.2 Program input 

The common units throughout the code are GeV and GeV/c^. The physical 
input parameters of GENTLE/4f an are defined in subroutine WUFLAG: 

ALFAI = a(0) = 137.0359895, the inverse running fine structure 

constant at zero momentum transfer 

ALFAS = as(lowg2) = 0.2 

ALPHS = as{2Mw) = 0.12 
ALPW = a{2Mw) = 1/128.07, the running fine structure con- 
stant at 2Mw 

AME = me = 0.51099906 X 10^3 GeV, the electron mass 

AMZ = Mz = 91.1888 GeV, the Z mass 

AMW = Mw = 80.230 GeV, the W mass 
GAMZ = Tz = 2.4974 GeV, the Z width 
GFER = Gf, = 1.16639 X 10^^ GeV-2, the Fermi constant □ 

For the fermion-boson couplings, the value a{2Mw) is used, while for the ISR 
corrections, the Thompson limit a(0) of the fine structure constant is applied. 
Depending on the flag value IINPT, the weak mixing angle is computed as 

SIN2W = sin^ 9w 
SIN2W = sin^ Ow 

SIN2W = sin^ Ow = 0.22591 (numerical user input) for IINPT = 2 (3.3) 

For charged current reactions, the following derived quantities are used in 

^ ag{2M\Y) is used for the QCD corrections to the W width and to the final 
state quark-anti-quark pairs in CC processes. The low momentum transfer value 
Q^sd^l 2Mvi/) is used for the NC32 process [19] in the calculation of gluon exchange 
diagrams. 



J- , ^ for IINPT = 3.1 

V2MIG, 

1 - M^r/Ml for IINPT = 1 (3.2) 
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GENTLE/4fan: 

GAMW = Tw= G.Ml, [l + 2^sWd^ for IGAMW = (3.4) 



GAMW = Tw = 2.085 (numerical user input) for IGAMW = 1 



(3.5) 



as well as 



GAE 



AswCw 



47ra(2Mi 



w 



AswCw 



GVE 



GAE( 1 - 4s 



w , 



GWF 



2V2 



-GAEa/2cw 



GWWG 



A'Ka{2Mw) 



GWWZ 



GWWG 



Sw 



GVE and GAE are the electron vector and axial vector couplings, GWF is the 
fermion-iy coupling, and GWWG and GWWZ are the triple gauge boson couplings 
for the photon and the Z boson in the Standard Model, respectively. All 
other couplings are uniquely determined by the above. Although, for histori- 
cal reasons, some of the variable names are slightly different, the same derived 
quantities are used for neutral current reactions. For the neutral current pro- 
cesses, the final state fermion types and their masses are part of the input. 
For the NCOS family of processes, this input is set in the subroutine NCINOO. 
The fermion types for the NC32 family are set by the flags I FERMI and IFERM2, 
while the corresponding masses are loaded from a DATA statement in the sub- 
routine bbmminQ. In the subroutine WUFLAG, the Higgs mass AMHIG and the 
Higgs width GAMHIG as well as the flags which steer the program flow are 
initialized. The Higgs contribution may be switched on or off by setting flag 
IHIGGS. 



The flags and, for the CC branch, the physical input may be overridden by 



Note that the fermion masses are stored in the variables AMI and AM2 for the NCOS 
as well as the NC32 families. This is, however, unproblematic, because the routines 
for each family are completely separate and are never executed together in a single 
run. 
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calls to the subroutine WUFLAG; see also section 3.5. We strongly recommend 
the user to carefully control that only options are active which are definitely 
declared to be compatible with each other. 
In WUFLAG, the following flags are initiahzed: 



I AND 
lANO: 

IANO=0: 

IAN0=±1 

IAN0=±2 

IAN0=±3 

IAN0=±4 

IAN0=±5 

IAN0=±6 

IAN0=±7 



IAN0=8: 



IBCKGR 
IBCKGR: 
IBCKGR=0: 

IBCKGR=1: 



to set the anomalous couplings defined in equations (B.3) 

to (B.8) of appendix B 

all couplings at their SM value 

XG = = ±0.5 , all other couplings at their SM values 
YG = 1/^ = ±0.5 , all other couplings at their SM values 
XZ = = ±0.5 , all other couplings at their SM values 
YZ = = ±0.5 , all other couplings at their SM values 
ZZ = zz = ±0.1 , all other couplings at their SM values 
DZ = 6z = ±0.5 , all other couplings at their SM values 
"HISZ scenario" [20] with 

X-y = ±0.1 , 

xz — Xj (cos^ 6w — sin^ 6w) / (2 sinOw cos^^^) , 
Sz — Xj/{2 sinOw cos9w) , 

and = yz = Zz = 

user-set values of x^, y^, Xz, yz, ^z, ^7, and Sz', 
for extra Z' couplings, see appendix C 



only active for IPR0C=1 

calculations for the double resonant CC03 process only, no 
ceil background 

inclusion of CCll background in accordance with ICHNNL 
not available for IDCS=1 



IBIN 
IBIN: 
IBIN=0: 
IBIN=1: 



IBINNU 
IBINNU: 
IBINNU=n: 



IBORNF 

IBORNF=0: 
IB0RNF=1: 



only active for IDCS=1 

calculation of dcr/dcos^^ for IBINNU fixed values of cos^ 
calculation of dcr/dcos6' integrated over each of the 
IBINNU cos6l-bins 



only active for IDCS=1 

number of points or bins for the differential cross-section 
calculation 



Born approximation, no radiative corrections 
initial state QED corrections are included 
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ICHNNL 



ICHNNL: 
ICHNNL=0: 

ICHNNL=1: 
ICHNNL=2: 

ICHNNL=3: 

ICHNNL=4: 
ICHNNL=5: 

ICOLMB 
ICOLMB: 



ICOLMB=0: 
IC0LMB=1: 

IC0LMB=2: 

IC0LMB=3: 
IC0LMB=4: 
IC0LMB=5: 

ICONVL 
ICONVL: 

ICONVL=0: 
IC0NVL=1: 



IPCS 

IDCS=0: 

IDCS=1: 



only active for IPR0C=1 and IBCKGR=1 

no background taken into account, calculation of inclusive 

CC03 quantities 

calculation of CC09 quantities for a leptonic final state 

calculation of CCIO quantities for a semi-leptonic final 

state with positively charged lepton 

calculation of CCIO quantities for a semi-leptonic final 

state with negatively charged lepton 

calculation of CCll quantities for a hadronic final state 

calculation of inclusive quantities for the CCll family 



determines, how the Coulomb correction is included in 
the differential cross-sections of equations (2.4) or (2.6) 
only relevant for the CC03 process 
Coulomb correction not included 

Coulomb correction as in equation (5) of reference [11], 
but without the phase factor 

Coulomb correction as in equation (5) of reference [11], 
but with A = (si — S2)/s' = 

Coulomb correction as in equation (5) of reference [11] 
Coulomb correction as in reference [14] 
Coulomb correction as in reference [12] 

only active for IPR0C=1 and IPR0C=2 

flux function convolution as given in equation (2.5) 
structure function convolution as indicated in eq. (2.10) 
attention: extremely CPU time consuming 

only active for IPR0C=1 and IBCKGR=0 

calculation of total cross-sections and moments 

calculation of the differential cross-section d(T/dcos6' with 

anomalous couphngs as determined by the fiag I AND 

if used in conjunction with ISR corrections, the structure 

function approach IC0NVL=1 must be used 

needs IBCKGR=0 
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IFERMl , IFERM2 
IFERMl: 



IFERM2: 



IGAMW 

IGAMW=0: 

IGAMW=1: 

IGAMZS 

IGAMZS=0: 

IGAMZS=1: 

IHIGGS 
IHIGGS: 
IHIGGS=0: 
IHIGGS=1: 

IINPT 
IINPT=0 
IINPT=1 
IINPT=2 

IiqCD 
IIQCD: 
IIQCD=0: 
IIQCD=1: 



IMOMN 
IMOMN: 
IMQMN=n: 



only active for IPR0C=2 IFERMl is an integer specifying 

the first final state fermion pair according to 

the PDG particle numbering scheme [13] 

like IFERMl, but for the second final state fermion pair 

remember that it is the second final state fermion pair 

that couples to the Higgs boson, if a calculation with the 

Higgs is requested 

computation of the W width according to (3.4) 
numerical user input 

constant Z width for Z propagators 
s-dependent Z width for Z propagators 



only active for IPR0C=2 
Higgs not included 
Higgs included 

weak mixing angle defined by equation (3.1) 
weak mixing angle defined by equation (3.3) 
numerical user input, only active for IPR0C=1 

only active for IPR0C=1 

no QCD corrections, i.e. ag{2Mw) = 

inclusive QCD corrections to the W width and final state 

QCD corrections to the final state quark pair are taken 

into account with ag{2MiY) — 0.12 



only active for IPR0C=1 and IDCS=0 

the n*^ moment of the physical quantity determined by 
the current I REG I M is computed 



The moments to be computed are determined by the values of IMMIN and 
IMMAX which control the scope of the Fortran loop DO IMOMN=IMMIN , IMMAX . 
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INCPRC 

INCPRC: 

INCPRC=0 

INCPRC=1 

INCPRC=2 

lONSHL 
lONSHL: 
IONSHL=0: 
I0NSHL=1: 

IPROC 

IPR0C=1 
IPR0C=2 
IPR0C=3 

IQEDHS 
IQEDHS: 



IQEDHS=0: 



IQEDHS=1 
IQEDHS=2 
IQEDHS=3 

IREGIM 
IREGIM: 



IREGIM=0: 

IREGIM=1: 
IREGIM=2: 

IREGIM=3: 

IREGIM=4: 



only active for IPR0C=3 

NC02 cross-section calculation, i.e. Z pair production only 
photon pair production only 
NCOS cross-section calculation 



only active for IPR0C=1 

on-shell limit for W pair production 

off-shell W bosons 



ceil family of processes 
NC32 family of processes 
NCOS family of processes with complete ISR 

determines the soft+virtual and hard radiators S and H to 
be used for universal corrections to total cross-sections 
only active, if IB0RNF=1 and ICONVL=0 
0{ol) exponentiated universal ISR corrections as given 
in equation (2.6) with the radiators S and H from 
equations (2.8) and (2.9) 

as IQEDHS=0, but terms of O {[oiLf] added to_S' and_ff 
as IQEDHS=1, but terms of O \c?L) added to_^ and_^ 
as IQEDHS=2, but terms of O (a^) added to S and E 

only active for IPR0C=1 and IDCS=0 

determines the physical quantity to be computed in addition 
to the total cross-section (see sections 2.3 and 2.4) 
computation of the total cross-section, 
used as normalization for next cases 

computation of moments of the radiative invariant mass loss 

computation of moments of the radiative energy loss 

for the structure function approach, i.e. for IC0NVL=1 only 

computation of moments of the W mass shift for the CC03 

process 

computation of W production angular moments for the CC03 
process with Chebyshev polynomials (compare section 2.4) 



The physical quantities to be computed are determined by the integer values 
of the variables IRMAX and IRSTP which control the scope of IREGIM via the 
Fortran statement DO IREGIM=0, IRMAX, IRSTP. 
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ITNONU 
ITNONU: 
ITNONU=0: 

ITN0NU=1: 



IZERD 
I ZERO: 
IZERO=0: 

IZER0=1: 

IZETTA 
IZETTA: 

IZETTA=0: 
IZETTA=1: 



only active for IB0RNF=1 

only universal hard ISR corrections are taken into 
account (compare equations (2.7) and (2.9) ) 
complete universal and non-universal hard ISR cor- 
rections are taken into account; only in effect for 
CC03 and for the NCOS family, IBCKGR=IDCS = 
attention: very CPU time consuming 



only for the flux function convolution, i.e. for ICONVL=0 
the constant term (Q;/7r)(7r^/3 — 1/2) in the S is neglected 
(compare equation (2.8) ) 

the soft-|-virtual S radiator is given by equation (2.8) 



only for the structure function convolution, i.e. for IC0NVL=1 

compare equation (4.21) of reference [5] 

ETA choice, i.e. in the hard photon radiation part of the 

structure function {2a/7T)L is used, whereas {2a/'n){L — 1) is 

entered into the soft part and the exponentiation 

BETA choice, i.e. {2a/7r){L — 1) is used for the whole structure 

function 



Finally, we mention that the relative and absolute precision requirements for 
the numerical integrations are set at different places in GENTLE/4f an. Sensi- 
ble default values for the precision requirements are chosen for the different 
branches of the program. In addition, we collect here the physical input pa- 
rameters which can be chosen by calls to subroutine WUFLAG: 



ALPHS, ALPW, AME, AMFER2, AMHIG, AMW, AMZ, DZ, GAMHIG, GAMW, GAMZ, GFER, 
SINW2, XG, XZ, YG, YZ, ZZ. 



3. 3 Routines of GENTLE/4f an 



In this subsection we give a brief description of the subroutines and functions 
in the code. GENTLE/4f an uses REAL*8 and C0MPLEX*16 variables. 

• FUNCTION alamk - computes the kinematic A-function A(x, y, z) — x'^ + y'^ + 

— 2xy — 2xz — 2yz 

• FUNCTION ALRUN - may be used in subroutine prop to compute photon 
propagators with running fine structure constant a 
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• SUBROUTINE bbmmin - constitutes the central routine for the 4f an part of 
the code, i.e. for computations for the NC32 family; performs the initializa- 
tion of physical and coupling constants for the 4f an part 

• FUNCTION c222 (c322, c322m, c422, c233) - computes the coupling 
function C222 (C322, C3225 ^422, C233) introduced in references [8,9]Q 

• FUNCTION cg422 (cg233) - computes the coupling function similar to C422 
(C233), but with a gluon exchange replacing the electroweak boson exchange 
in the final state (compare figure le of reference [19]) 

• SUBROUTINE check - checks and prints the input to the 4f an part; produces 
the integer error fiag ierr; prints messages if input errors are encountered 

• FUNCTION COULOM - computes, according to the fiag value of ICOLMB, the 
term COULOM for the Coulomb correction factor (1+COULOM) 

• FUNCTIONS DDILOG, DLI2 - provide the dilogarithm function Li2(x) for real 
argument; if x > 1, x is understood to be the real part of the argument 
with infinitesimal imaginary part, and the real part of the function value is 
returned 

• SUBROUTINE DELTPR - computes, for the fiux function approach, the three- 
fold (with respect to Si, S2, and s') differential, non-universal contributions 
to the CC03 cross-section 

• SUBROUTINE DELVNU - computes the non-universal soft+virtual contribu- 
tion to the NCqed cross-section 

• FUNCTION DLI3 - provides the trilogarithm function Li3(a;) for x < 1 

• FUNCTION DS12 - provides the Nielsen function Si,2(a;) for a; < 1 

• SUBROUTINE dsdshsz - provides the twofold differential cross-section 
d^cr/( dsi ds2) for the NC32 family as given in references [8,9] 

• FUNCTION DXS2 - provides, for the CC/NC branch of the code and according 
to the flags lONSHL, IPROC, IBORNF, IBCKGR, and ICHNNL, twofold differential 
expressions (with respect to Si and S2) for the cross-section or for moments 
by calling XSECB and XSECBG or XSEC2E 

• FUNCTION FACINC - is the inverse function of FACNC 

• FUNCTION FACINV - represents the inverse function of FACT 

• FUNCTION FACNC - is the mapping function for the s' integrations used in 
calls to FDSIMT 

• FUNCTION FACT - is the mapping function for the s', x^, and x_ integrations 
used in calls to FDSIMP and FDSIMS 

• FUNCTION FCONC - provides the integrand of the s' integration in XSNC2E 

• FUNCTION FCONVl - provides, according to the flag ICONVL, threefold dif- 
ferential (with respect to si, S2, and s' or x+) expressions for 

(1) the ISR-corrected cross-section, the radiative invariant mass loss, 
and the W mass shift using the flux function approach 

(2) total and differential cross-sections as well as moments in the struc- 
ture function approach by calling FC0NV2 



Beware: The normalization differs from the references in some cases. 
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• FUNCTION FC0NV2 - provides fourfold differential (with respect to Si, S2, 
x^, and a;_) expressions for all quantities available in the structure function 
approach 

• SUBROUTINE FDSIMP - is similar to SIMPS, but performs the integration 
with the help a mapping function which is called in the argument list 

• SUBROUTINES FEYNINTl, FEYNINT2, FEYNINT3 - compute, from the ana- 
lytical formulae of reference [6], the values of several complicated integrals 
needed for the non-universal hard ISR correction 

• FUNCTION FINE12 (FI12NC) - provides the differential cross-section 
da/ dxi, Xi — Si/s for the CC/NC (NCqed) branch of the code 

• FUNCTION FINE34 (FI34NC) - provides the differential cross-section 
do"/( dxi dx2) for the CC/NC (NCqed) branch of the code by calling the func- 
tion DXS2 (XSBRNC and XSNC2E) which actually provides the value 

• FUNCTION FMAPSH - performs, for the NC branch, the mapping in the invari- 
ant mass squared S2 of the fermion pair defined by the flag IFERM2 

• FUNCTION FMAPSZ - performs, for the NC branch, the mapping in the invari- 
ant mass squared Si of the fermion pair defined by the flag IFERMl 

• FUNCTION FMBW - maps Breit-Wigner resonance peaks 

• FUNCTION FMQED - maps peaks of the type l/s 

• FUNCTION FOTF - computes the leptonic vacuum polarization for ALRUN from 
lepton and effective quark masses 

• FUNCTION FOTFl - computes the leptonic vacuum polarization for ALRUN 
from lepton masses and the hadronic vacuum polarization of reference [21] 

• FUNCTION FSZ - performs, for the NC branch, the integration over S2; this 
is the complement to the function FINE 12 

• FUNCTION FTBINV - represents the inverse function of FXTB 

• FUNCTION FXSECB - provides, in dependence on the flag IDCS and for the 
CC03 process, the threefold differential cross-section d'^o"/( dxi dx2 d cos^^) 
or the argument T„(cos6';a6) d^(Tcco3/( dsi ds2 dcos6') of the moment inte- 
gral (2.15) 

• FUNCTION FXSECD - provides, in dependence on the flags IDCS and IBIN and 
for the CC03 process, the integral over the threefold differential cross-section 
d^(T/( dxi dx2 d cos^) for given limits 

• FUNCTION FXTB - is the mapping function for the cos^ integration with the 
subroutine FDSIMV 

• FUNCTION g222 (g322, g322m, g422, g233) - computes the kinematic 
function ^222 (^322, ^322) ^422, fes) introduced in references [8,9] 

• SUBROUTINE GFBGD2 - computes the kinematic functions Qk for the pure 
background for the CCll process 

• SUBROUTINE GFBGD3 - computes the kinematic functions Qk for the signal- 
background interferences of the CCll process 

• SUBROUTINE GFBRN - computes the kinematic functions for the CC03 
process (compare equation (2.4), see reference [7] ) 

• SUBROUTINE GFBRNA - computes, in dependence on the flags IDCS and lANO, 
the angular kinematic functions ^^'^°^(<s; Si, S2, cos6') (compare eq. (2.13), see 



21 



reference [5] ) 

• SUBROUTINE GFBRNB - computes, in dependence on the flags IDCS and lANO, 
the integrated angular kinematic functions ^^"^"^(s; Si, S2, 0086*111111, cos6'iiiax) 

• SUBROUTINE GFBRNC - computes the kinematic function for the NCqed branch 

• MAIN - performs the ramification between the NCqed and the CC/NC branches 
of the code; carries out some initiahzation; calls WWINOO and WWINOl for 
further initialization; loops over the center of mass energies; loops over the 
desired values of IREGIM and IMOMN; prints results 

• FUNCTION MYEXPl (MYEXP2) - is the inverse function of MYLOGl (MYL0G2) 
as physics input 

• FUNCTION MYLOGl (MYL0G2) - is the mapping function for the Si(s2) inte- 
gration used in calls to FDSIMP (FDSIMS) 

• SUBROUTINE NCGENT - the central routine for the NCqed branch of the code; 
sets integration precisions; calls NCINOO and NCINOl for further initiahzation; 
loops over the center of mass energies; provides the total cross-section for 
the NCqed branch; prints results 

• SUBROUTINE NCINOO - initializes the flags needed for the NCqed branch of 
the code; initializes the fermion types and couplings and performs other 
kinematics independent initialization for the NCqed branch 

• SUBROUTINE NCINOl - performs, for the NCqed branch of the code, the ini- 
tialization of quantities that depend on s only and provides the soft-|-virtual 
radiator 

• SUBROUTINE NCIN02 - provides the phase space factor and the coupling 
function for the NCqed branch given in equation (4.5) of first of references [6] 

• SUBROUTINE prop - provides the Z, 7, and gluon propagators for use in the 
subroutine dsdshsz 

• SUBROUTINE PROPAG - computes the neutral gauge boson propagators used 
in NCIN02 

• FUNCTION RHOINV - is the inverse function of RHOPR 

• FUNCTION RHOPR - provides, for the CC/NC branch of the code and according 
to the mapping choice via IMAPPG, the mapping function for the Si and S2 

integrations 

• FUNCTION RI3 ~ computes the scalar two-point function needed in FOTF and 
FOTFl 

• SUBROUTINE SETANO - initiahzes the anomalous couplings defined by equa- 
tions (B.3) to (B.8) in dependence on the fiags IDCS and lANO 

• SUBROUTINE SHORTI - computes the kinematic functions ^ccii'^ ^-nd ^ccii 
as given in equations (3.9) and (3.10) of reference [5] 

• SUBROUTINE SIMPS, SIMPT - performs the integration of an input function 
inside given limits and with a given precision by applying a self-adaptive 
Simpson algorithm 

• FUNCTION SPENC - is an auxiliary function used by the complex dilogarithm 
function XSPENZ 

• FUNCTION STRUCF - computes, according to the value of IZETTA, the struc- 
ture function D{x) from appendix A of reference [15] 
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FUNCTION TOTXS - computes the total cross-section for the CC/NC branch 

and the differential cross-section for the CC03 case 

FUNCTION TRILOG - provides the trilogarithm function Li3(a;) for < a; < 1 
FUNCTION TRIS12 - provides the Nielsen function Si,2(a;) for < a; < 1 
SUBROUTINE WUFLAG - initiahzes and changes the default flags as well as 

physics inpTit 

SUBROUTINE WWINOO - performs the kinematics independent initialization 
of mathematical and physical constants for the CC/NC branch of the code; 
calls the subroutines bbmmin and SETANO for further initialization 
SUBROUTINE WWINOl - performs, for the CC/NC branch, the initialization 
of quantities that depend on the center of mass energy or, equivalently, s 
only; performs initialization in dependence on the flags IPROC, ICONVL, and 
IQEDHS; provides the soft-|-virtual radiator for the CC/NC computation 
SUBROUTINE WWIN02 - performs the initialization of tree level quantities 
that depend on s, si, and S2; for the CC/NC branch 

SUBROUTINE WWIN03 - performs the initialization of radiative quantities that 

depend on s, Si, and S2; for the CC/NC branch 

SUBROUTINE XSBRNC - provides the tree level cross-section for the NCqed 
branch of the code 

SUBROUTINE XSEC2E - is needed for the computation of ISR-corrected cross- 
sections, moments, and distributions; performs the s' integration for the 
flux function approach and the x+ integration for the structure function 

approach 

SUBROUTINE XSEC3B - computes, in the flux function approach and for the 
ecu family with IBCKGR=1, the universally ISR-corrected, threefold (with 
respect to si, S2, and s') differential cross-section for background contribu- 
tions by calling the subroutine XSECBG 

SUBROUTINE XSEC3E - is used in the flux function convolution for the 
CC03 process, i.e. for IPR0C=1, IB0RNF=1, and ICONVL=0 only; computes, 
in dependence on the settings of ITNONU and IQEDHS and according to for- 
mula (2.6) and reference [10], the threefold differential CC03-expressions 

SUBROUTINE XSECB - for the CC/NC branch; provides, in dependence on the 
flags IPROC, IOCS, and IREGIM, the twofold differential (with respect to Si 
and S2) cross-section or moments for the chosen tree level process by calling 
GFBRN or integrating FXSECB over the W production angle 
SUBROUTINE XSECBG - provides, in dependence on the channel chosen via 
the flag ICHNNL, the twofold (with respect to Si and S2) differential cross- 
section for the signal-background interference and the pure background in 
the ecu process by calling the subroutines GFBGD2 and GFBGD3 
SUBROUTINE XSNC2E - provides the twofold (with respect to si and S2) dif- 
ferential cross-section for the radiative NCqed cross-section by integrating 
over s' 
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• SUBROUTINE XSNC3E - computes, according to the flags ITNONU, the three- 
fold (with respect to Si and S2) differential, radiative NCqed cross-section 

• FUNCTION XSPENZ - provides the complex dilogarithm function Li2{z) by 
transforming the argument to the area of fast convergence 

3.4 Program output 

The output of GENTLE/4fan strongly depends on the user-chosen branch of 
the program or, in other words, on the flag settings. Four different output 
appearances are to be distinguished, namely 

(1) the output of total cross-sections and moments for the CCll family 

(2) the output of angular distributions for the CC03 process 

(3) the output of the 4f an branch of the code 

(4) the output of the NCqed branch of the code 

In the below subsections, we will provide and explain examples for each of the 
above four kinds of output. The below outputs will, at the same time, serve 
as test run output, where we assume that the user has used the flags and 
physical input as shown in the output. Any other pre-set input should remain 
unchanged for test run purposes. 

3.4-1 Output of total cross-sections and moments for CCll 

The output of GENTLE/4f an for the total cross-section and moments for the 
CCll family of processes appears as follows. 

***** This is GENTLE/4fan — Version 2.0 ***** 

This is the CC branch of GENTLE/4fan for total cross-sections 
FLAGS :IPROC ,IINPT , lONSHL, IBORNF , IBCKGR, ICHNNL= 111112 
FLAGS : IGAMZS , IGAMW , ITNONU, IQEDHS,ICOLMB,ICONVL= 10 12 
FLAGS :IZERO , IZETTA, IIQCD , IOCS ,IANO ,IBIN =011000 
FLAGS :IMAP , IRMAX , IRSTP , IMMIN ,IMMAX =11114 
EPSl = l.OOOOOOOOOOOOOOOOE-04 
S1MIN,S1MAX= O.OOOOE+00 0.2592E+05 
S2MIN,S2MAX= O.OOOOE+00 0.2592E+05 
Energy (GeV) Cross-section (pb) 

161.00 0.1336661D+00 
MOMENTS 

11 2 3 4 

1 0.4753809D+00 . 3387719D+01 . 5607608D+02 0. 1411261D+04 
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O.OOOOE+00 0.3063E+05 
O.OOOOE+00 0.3063E+05 
Cross-section (pb) 
0.4952086D+00 



S1MIN,S1MAX= 
S2MIN,S2MAX= 
Energy (GeV) 

175.00 
MOMENTS 

11 2 3 4 

I 0.1126753D+01 . 8432371D+01 . 1030576D+03 0. 1885898D+04 
S1MIN,S1MAX= O.OOOOE+00 0.3610E+05 

S2MIN,S2MAX= O.OOOOE+00 0.3610E+05 
Energy (GeV) Cross-section (pb) 

190.00 0.6080011D+00 
MOMENTS 

II 2 3 4 

I 0.2151788D+01 . 2755849D+02 . 4972633D+03 0. 1089393D+05 
S1MIN,S1MAX= O.OOOOE+00 0.4203E+05 

S2MIN,S2MAX= O.OOOOE+00 0.4203E+05 
Energy (GeV) Cross-section (pb) 

205.00 0.6355688D+00 
MOMENTS 

II 2 3 4 



0.3208536D+01 . 5908307D+02 . 1485244D+04 . 4348380D+05 



***** This is GENTLE/4fan — Version 2.0 ***** 



This is the CC branch of GENTLE/4fan for total cross-sections 



FLAGS :IPROC ,IINPT , lONSHL , IBORNF , IBCKGR, ICHNNL= 
FLAGS : IGAMZS , IGAMW , ITNONU , IQEDHS , ICOLMB , ICONVL= 
FLAGS :IZERO , IZETTA, IIQCD ,IDCS ,IANO ,IBIN = 
FLAGS :IMAP ,IRMAX ,IRSTP ,IMMIN ,IMMAX 
EPSl = l.OOOOOOOOOOOOOOOOE-04 
S1MIN,S1MAX= O.OOOOE+OO 0.2592E+05 
S2MIN,S2MAX= O.OOOOE+OO 0.2592E+05 
Energy (GeV) Cross-section (pb) 

161.00 . 1332829D+00 

MOMENTS 

3 1 2 3 

3 -0.3052332D+01 . 2605644D+02 -0 . 3833234D+03 
S1MIN,S1MAX= O.OOOOE+00 0.3063E+05 
O.OOOOE+00 0.3063E+05 
Cross-section (pb) 
0.4946375D+00 



S2MIN,S2MAX= 

Energy (GeV) 
175.00 

MOMENTS 
3 1 
3 -0.5666559D+00 



2 3 
0.8922053D+01 -0. 95531 19D+02 



1 1 



1 1 
3 3 



10 2 

12 



1 4 



0.7596405D+04 



0.1918868D+04 
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S1MIN,S1MAX= O.OOOOE+00 0.3610E+05 

S2MIN,S2MAX= O.OOOOE+00 0.3610E+05 

Energy (GeV) Cross-section (pb) 

190.00 0.6074997D+00 
MOMENTS 

3 1 2 3 4 

3 -0.1004648D+00 . 9783632D+01 -0 .4987223D+02 0. 1646710D+04 

S1MIN,S1MAX= O.OOOOE+00 0.4203E+05 

S2MIN,S2MAX= O.OOOOE+00 0.4203E+05 

Energy (GeV) Cross-section (pb) 
205.00 0.6352834D+00 

MOMENTS 

3 1 2 3 4 

3 0.1802810D+00 . 1232392D+02 . 6206403D+01 . 2242722D+04 



First, the program identifies itself and tlic brancfi of the code for which output 
is produced. Then, all relevant flags are printed, before the relative precision 
EPSl of the last computation is given. Finally, physics output is generated. 
For each center of mass energy the energy and the corresponding total cross- 
section are printed. Then, in the first column, integers k indicate the value 
of IREGIM and thus determine the physics quantity for which moments are 
output. The first line in the moments' output represents the degree n of the 
moments in the column below. In the second line, the moment of degree n 
for the quantity determined by k is given. Its value is normalized to the total 
cross-section. 



3.4-2 Output of differential cross-sections for CC03 

First we give an example of output for IBIN=0, with differential cross-sections 
computed at fixed values for the scattering angle. 

***** This is GENTLE/4fan — Version 2.0 ***** 

This is the CC branch of GENTLE/4fan for diff . cross-sections 

XG= O.lOOOOD+00 XZ= -O.lOOOOD+00 

YG= O.OOOOOD+00 YZ= O.OOOOOD+00 

ZZ= O.OOOOOD+00 DZ= 0.50000D-01 
FLAGS :IPROC ,IINPT , lONSHL, IBORNF, IBCKGR, ICHNNL= 1110 2 
FLAGS :IGAMZS,IGAMW , ITNONU, IQEDHS, ICGLMB, ICONVL= 10 12 1 
FLAGS :IZERO , IZETTA, IIQCD , IDCS ,IANO ,IBIN =011180 
FLAGS :IMAP ,IRMAX ,IRSTP ,IMMIN ,IMMAX =13 5 11 

EPSl = l.OOOOOOOOOOOOOOOOE-04 
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S1MIN,S1MAX= O.OOOOE+00 . 2592E+05 
S2MIN,S2MAX= O.OOOOE+00 . 2592E+05 
Energy (GeV)= 161.0000000000000 
COSW Cross-section (pb) 
-1.000 0.0490713 
-0.800 0.0536632 
-0.400 0.0646671 
0.000 0.0794040 
0.400 0.1009377 
0.800 0.1381398 
1.000 0.1760354 
S1MIN,S1MAX= O.OOOOE+00 . 3063E+05 
S2MIN,S2MAX= O.OOOOE+00 0.3063E+05 
Energy (GeV)= 175.0000000000000 
COSW Cross-section (pb) 
-1.000 0.0951087 
-0.800 0.1123347 
-0.400 0.1564827 
0.000 0.2238638 
0.400 0.3411499 
0.800 0.5873924 
1.000 0.8353949 
S1MIN,S1MAX= O.OOOOE+00 0.3610E+05 
S2MIN,S2MAX= O.OOOOE+00 0.3610E+05 
Energy (GeV)= 190.0000000000000 
COSW Cross-section (pb) 
-1.000 0.0710393 
-0.800 0.0891843 
-0.400 0.1359378 
0.000 0.2124469 
0.400 0.3669140 
0.800 0.7949695 
1.000 1.3546129 
S1MIN,S1MAX= O.OOOOE+00 0.4203E+05 
S2MIN,S2MAX= O.OOOOE+00 0.4203E+05 
Energy (GeV)= 205.0000000000000 
COSW Cross-section (pb) 
-1.000 0.0512161 
-0.800 0.0678664 
-0.400 0.1099483 
0.000 0.1806793 
0.400 0.3368272 
0.800 0.8655771 
1.000 1.7575708 
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After the identification of the active GENTLE/4f an branch, the anomalous cou- 
phngs (compare page 15 and appendix B) used in the CC03 subprocess are 
output, before the used flags and the relative precision of the run are printed. 
Next, below the print of each center of mass energy, the differential cross-sec- 
tions for IBINNU fixed values of the scattering angle COSW are written to the 
output. 

Below, we present output for IBIN=1, i.e the given differential cross-section 
values are integrated over bins in the scattering angle. The appearance of the 
output is identical to the above output for IBIN=0 except that, for each energy, 
the total cross-section is given right below the prints of scattering angle values 
and the corresponding cross-sections. 

***** This is GENTLE/4fan — Version 2.0 ***** 

This is the CC branch of GENTLE/4fan for diff . cross-sections 

XG= O.lOOOOD+00 XZ= 0.65543D-01 

YG= O.OOOOOD+00 YZ= O.OOOOOD+00 

ZZ= O.OOOOOD+00 DZ= 0.11957D+00 
FLAGS :IPROC ,IINPT , lONSHL, IBORNF, IBCKGR, ICHNNL= 1110 2 
FLAGS:IGAMZS,IGAMW , ITNONU, IQEDHS , ICOLMB , ICONVL= 10 12 1 
FLAGS :IZERO , IZETTA, IIQCD , IDCS ,IANO ,IBIN =011171 
FLAGS :IMAP ,IRMAX ,IRSTP ,IMMIN ,IMMAX =13 5 11 

EPSl = 1 . OOOOOOOOOOOOOOOOE-04 
S1MIN,S1MAX= O.OOOOE+00 0.2592E+05 
S2MIN,S2MAX= O.OOOOE+00 0.2592E+05 
Energy (GeV)= 161.0000000000000 

COSW Cross-section (pb) 

-1.000 - -0.600 0.020923 
-0.600 - -0.200 0.025547 
-0.200 - 0.200 0.031763 

0.200 - 0.600 0.040919 

0.600 - 1.000 0.057338 
TOTAL XSEC = 0.176489 
S1MIN,S1MAX= O.OOOOE+00 0.3063E+05 
S2MIN,S2MAX= O.OOOOE+00 0.3063E+05 
Energy (GeV)= 175.0000000000000 

COSW Cross-section (pb) 

-1.000 - -0.600 0.042294 
-0.600 - -0.200 0.060777 
-0.200 - 0.200 0.089224 

0.200 - 0.600 0.139347 

0.600 - 1.000 0.247555 
TOTAL XSEC = 0.579198 
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S1MIN,S1MAX= 0, 
S2MIN,S2MAX= 0. 
Energy (GeV)= 

COSW 
-1.000 - -0.600 
-0.600 - -0.200 
-0.200 - 0.200 
0.200 - 0.600 
0.600 - 1.000 
TOTAL XSEC = 
S1MIN,S1MAX= 0. 
S2MIN,S2MAX= 0. 
Energy (GeV) = 

COSW 
-1.000 - -0.600 
-0.600 - -0.200 
-0.200 - 0.200 
0.200 - 0.600 
0.600 - 1.000 
TOTAL XSEC = 



OOOOE+00 0.3610E+05 
OOOOE+00 0.3610E+05 
190.0000000000000 
Cross-section (pb) 

. 032453 

0.051912 

0.084241 

0.150952 

. 346228 

. 665786 
OOOOE+00 0.4203E+05 
OOOOE+00 0.4203E+05 
205.0000000000000 
Cross-section (pb) 

. 023983 

0.041371 

0.071246 

0.139243 

0.393478 

. 669320 



3.4-3 Output for the 4:13.11 branch 
***** This is GENTLE/4fan — Version 2.0 ***** 



This is the 4fan branch of GENTLE/4fan 
4fan was called for fermions 13 5 
with masses 0.500000D-03 O.lOOOOOD-01 
total cross section is calculated 



FLAGS :IPROC , IINPT , lONSHL, IBORNF , IBCKGR, ICHNNL= 2 10 
FLAGS : IGAMZS , IGAMW ,ITNONU,IQEDHS,ICOLMB,ICONVL= 10 10 
FLAGS :IZERO ,IZETTA,IIQCD , IOCS ,IANO ,IBIN =110 
FLAGS :IMAP ,IRMAX , IRSTP , IMMIN , IMMAX =10 111 

FLAGS FOR 4FAN: IFERMl, IFERM2, IHIGGS=13 5 
EPSl = l.OOOOOOOOOOOOOOOOE-04 
ABSl = 1 . OOOOOOOOOOOOOOOOE-04 
S1MIN,S1MAX= 0.5805E+04 0.1128E+05 
S2MIN,S2MAX= 0.9000E+03 0.2592E+05 
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Energy (GeV) 

161.00 
S1MIN,S1MAX= 
S2MIN,S2MAX= 
Energy (GeV) 

175.00 
S1MIN,S1MAX= 
S2MIN,S2MAX= 
Energy (GeV) 

190.00 
S1MIN,S1MAX= 
S2MIN,S2MAX= 
Energy (GeV) 

205 . 00 



Cross-section (pb) 

0.2907564D+00 
0.5805E+04 0.1128E+05 
0.9000E+03 0.3063E+05 
Cross-section (pb) 

0.8156756D+00 
0.5805E+04 0.1128E+05 
0.9000E+03 0.3610E+05 
Cross-section (pb) 

0.9510082D+01 
0.5805E+04 0.1128E+05 
0.9000E+03 0.4203E+05 
Cross-section (pb) 

. 1249948D+02 



3.4-4 Output for the NCqed branch 

***** This is GENTLE/4fan — Version 2.0 ***** 

This is the NCqed branch of GENTLE/4fan 
FLAGS : NINPT , NNCPRC , NGAMZS , NBORNF , NTNONU , NQEDHS 
2 110 

GFER = 1 . 166390000000000E-05 
AMW = 80.23 
AMZ = 91.1888 
GAMZ = 2.4974 

ALPW = 7 . 808229874287499E-03 
SIN~2(\theta_W) = .2310309124515784 

IFERMl , IFERM2 =2 4 
RNC0U1,RNC0U2 = 1.00000000 3.00000000 
AMI ,AM2 = .105658389 4.3 

CUTM12,CUTM34 = .211316778 8.6 
CUTXPR = .0 

EPSl l.OOOOOOOOOOOOOOOE-05 
ABSl l.OOOOOOOOOOOOOOOE-15 
ENERGY (GeV) CROSS-SECTION (fb) 

161.0 .8105861091D+02 

175.0 .6725488126D+02 

190.0 .6303702520D+02 

205 . . 5842607173D+02 
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Table 3 

Correspondence between fermion pair indices and fermion types for the NCqed branch 
o/GENTLE/4fan. 



F(;rmi(ni Pair Index 


1 


2 


3 


i 


Fermion Type 






u, c 


d, s, b 



After identifying the active branch, the flags relevant to the NCqed branch of 
the code arc printed. The flags are as set in WUFLAG. but the flrst letter is 
changed to an N standing for NCqed. Next, the relevant physical input con- 
stants are printed the meaning of which is described in section 3.2. The next 
paragraph of output begins with the fermion pair indices. The correspondence 
of fermions and indices is given in table 3. Then, for the two fermion types used 
in the computation, the color factors (RNCOUl, RNC0U2), the fermion masses 
(AMI, AM2), and the lower invariant mass cuts (CUTM12, CUTM34) are printed. 
Closing this paragraph of output, the minimum value CUTXPR of the fraction 
s'/s is printed. Finally, together with the relative and absolute required pre- 
cisions EPSl and ABSl, the total cross-sections are output for each required 
center of mass energy. 

3.5 Use of the program 

In this section, we will give a short guide to the compilation and use of 
GENTLE/4fan. 

We recommend to exhaust all means of compiler optimization, because, ac- 
cording to our experience, the CPU time consumption can thus be reduced 
by a factor of two or more. When available, the use of a compiler option that 
makes the default size of floating-point constants REAL*8 may be advanta- 
geous. For HP computers, we recommend 

f77 +04 +Onolimit +ppu -K -R8 gentle. f . 
For Silicon Graphics computers, compilation with 

f77 -r8 -03 -non_shared -G 32 -jmpopt -mips2 -static gentle. f 
has proved successful. 

Unless the user wants to run GENTLE/4f an with default flags, the user should 
set the flags described in section 3.2. The setting of a flag IFLAG to a value 
IVALUE is accomplished by adding the hne 
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CALL WUFLAG ( ' IFLAG ' , I VALUE) 



after the DATA statements in the MAIN of GENTLE/4f an. Please note that the 
numerical format of IVALUE is INTEGER. Similarly, for IAN0=8, the anomalous 
couphngs XG, YG, XZ, YZ, ZZ, DG, and DZ may be set by calls hke 

CALL WUFLAG ( ' XG ' , VALUE) 

at the same place in the MAIN. Please note that the numerical format of VALUE 
is DOUBLE PRECISION. In addition, the physical input variables GFER, ALPW, 
AME, AMW, AMZ, GAMZ, and ALPHS for the CC branch of the code as well as AMHIG 
and GAMHIG for the 4f an branch can be changed by analogous calls to WUFLAG. 
Finally, it is possible to change the physics input to the 4fan and NCqed 
branches of the code in the subroutines bbmmin and NCINOO respectively. 

To determine the physics quantities and the degrees of moments to be com- 
puted and printed, the user must set the range of the loops over IREGIM and 
IMOMN (see page 18). 

Cuts on the final state invariant pair masses may be applied in the CC and the 
4f an branches of the code by altering SIMIN, SIMAX, S2MIN, and S2MAX in the 
MAIN. In the NCqed branch, final state invariant pair mass cuts are implemented 
by setting X12MIN, X12MAX, X34MIN, and X34MAX in the subroutines NCGENT and 
FINC34 ( X12 = si/s, X34 Ss/s). 
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A Numerical Integration in GENTLE/4f an 



The integrations over si,S2 are performed numerically by a successive use 
of an adaptive one-dimensional Simpson integration routine. One (two) more 
integrations are needed to calculate initial state QED corrections in the fiux 
function (structure function) approach. Mappings of the integrands may be 
chosen with flag IMAP. We now describe the mappings. 

After the analytical integration over the six angular variables, one has to treat 
Breit-Wigner resonances arising from the exchange of intermediate bosons: 

In the case of neutral current processes, additional peaks appear at small 
momentum transfer due to virtual photons or gluons, 

fiow{V,s)^-, V^^,g. (A.2) 

s 

For charged current processes, the standard mapping of the Breit-Wigner res- 
onance of the W boson is applied in GENTLE/4f an. 

In the case of neutral current processes, the peaks at small momentum transfer 
and the Brcit-Wigncr resonances of the Z and the Higgs boson have to be 
mapped simultaneously. In GENTLE/4f an, the singularity at small momentum 
transfer is mapped first by the usual transformation to a new variable s, 

s = /i(s) = lns. (A.3) 

Then follows a mapping of the Z resonance, which does not destroy the pre- 
vious mapping, 

= C S - My . . ^. 

s — 12(8) — s + 7= — r-— r arctan -= — 7-— r- (A.4) 

The variables My and [FyMy] are calculated taking into account the previous 
mapping. 

Ml = A(M^), [VyMy] = h{Ml + VyMy) - MM^). (A.5) 

The function /2 is inverted numerically by Newton's Method, 

Si+i = - 7773T, (A.6) 
which always converges for the starting value sq — My. The free constant c 
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in equation (A. 4) is optimized for every final state separately. However, the 
quality of the mapping (A. 4) is sensitive to the order of magnitude of c only. 
In the case of Higgs production, a procedure similar to (A. 4) is added to map 
the Higgs resonance. 



B Anomalous Triple-Boson Couplings 



In addition to the SM Lagrangian, we consider the following CP conserving 
operators of dimension 6 [^: 

A£,ff = g'^ {D,^)^ B^^ {D,^) + g^ ri^^^ 



In the unitary gauge, these operators lead to the following effective Lagrangian 
for the WWV vertex: 



rWWV _ • 

eff ~ ^9wwv 



9i {w;^uW-^ - w-^^w-,) 

^ (B.2) 



where V G {Z, 7}. Therefore, one has six relevant anomalous triple-boson 
couphngs, namely 5(7) 9i 1 '^7; i^z, \, and \z- Electromagnetic gauge invariance 
requires gj = 1. In the HISZ-scenario [20] (IANO= ±7), one sets a^^i = c^vy^i- 
Thus, in this scenario, one obtains relations between the parameters k^, kz, 
and gf . 

For the WWZ-vertex, we also include a C and P violating term in the La- 
grangian: 

£z = - ^ daZp^ (w^" dP W-" - W^" dP iy-°) , 

where we have introduced the dual field strength tensor 
7 — —1= 7"^ 

To disentangle the anomalous couplings' contributions to the cross-section, we 
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For definitions of anomalous couplings, see e.g. [22,16]. 
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use redefined parameters: 



5z = {9f -I) cot Ow, (B.3) 

X-y = K-y — 1, (B.4) 

xz = (kz - 1) (cot ew + 5z), (B.5) 
^7 = (B-6) 
yz = Azcot^i^, (B.7) 

zz = zz. (B.8) 

The GENTLE/4f an user can set values for the anomalous couplings defined in 
equations (B.3) to (B.8). In the SM, the anomalous couplings (B.3) to (B.8) 
vanish. The special case of anomalous couplings induced by a heavy extra 
neutral gauge boson is described in appendix C. 

C Treatment of Extra Neutral Gauge Bosons 

The effect of an extra neutral gauge boson Z' can be described with GENTLE/- 
4f an by two anomalous couplings Qy^wi ^^'^ 9wwz of the photon and the 
Standard Model Z boson to W pairs [23], 

Qww-i = 1 + (^7 and g^wz = cot 9w + (^z, (C.l) 
where 

(^7 = 9wwzi ( — - — ) ^(1 + Ax)xz + gwwz2 ( — - — ) VX2, 
\ a V J \ a V J 

5z = -cotew + gwwzi — 0- + ^x)+gwwz2 ——, (C.2) 

a a Xz 

and 

s . _ 2Mz{Mz - Ml) s 

In equations (C.2) and (C.3), Mj, f j, Oj, i = 1,2 denote the masses and electron 
couplings of the mass eigenstates Zi and Z2 , which are, in general, a result 
of a mixing of the gauge symmetry eigenstates Z and Z', 
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Z\ — Z COS Qm + "Zi' sin ^m, -^2 = — sin Qu + -Z^' cos ^m, 
= V cos + "v' sin 6*^, '^'2 = —"v sin + f' cos du^ 
ai = acos^M + o'sin^M, 02 = — asin^M + o'cos^m, (C.4) 

where v, a (v', a') denote the couphngs of the symmetry eigenstates Z {Z') to 
electrons and 9m is the ZZ' mixing angle. The couplings gwwzi and gwwz2 
are fixed by the condition that only the Z couples to W pairs, 

9wwzi = cot 9w cos 9m, gwwz-z = cot 9w sin 9m- (C.5) 

The mass shift from Mz to Mi in the Zi propagator is absorbed into 6^ and 
6z through Ax assuming - Mi < Mz and - Mi < (s - M|) /(2Mz). 
Both approximations are equivalent for W pair production. 



D Options of GENTLE 

Table D.l shows the observables GENTLE can calculate and the options, which 
may be chosen. Please note that in the CCll family anomalous couphngs are 
included only for the CC03 process. 



Table D.l 

Options of GENTLE; SF - structure function approach; FF - flux function approach; 
h.o. - higher order QED corrections; Coul. - Coulomb singularity. 





IREGIM 


on-shell 


Born 


SF 


FF 


h.o. 


Coul. 


non-uni. 


anom.c. 


^CC03 





+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


^CC03 


1 


+ 




+ 


+ 


+ 


+ 


+ 


+ 


^CC03 


2 


+ 




+ 




+ 


+ 




+ 


i^ccos 


3 


+ 




+ 


+ 


+ 


+ 


+ 


+ 




4 


+ 


+ 


+ 




+ 


+ 






^^CC03 

d cos^ 




+ 


+ 


+ 




+ 


+ 




+ 


^CCll 







+ 


+ 


+ 


+ 


+ 




+ 


^CCll 


1 






+ 


+ 


+ 


+ 




+ 


^CCll 


2 






+ 




+ 


+ 




+ 


^CCll 


3 






+ 


+ 


+ 


+ 




+ 








+ 




+ 


+ 




+ 




I^NCOS 






+ 




+ 


+ 




+ 




iyNC32 






+ 


+ 


+ 


+ 









36 



References 



[1] D. Bardin, R. Kleiss et al., Report of the Working Group on Event Generators 
for WW Physics, in G. Altarelli, T. Sjostrand, F. Zwirner, eds., Proceedings of 
the Workshop on Physics at LEP 2, CERN Yellow Report CERN 96-01 (1996), 
Vol. 2, 3-101, and references therein. 

[2] F. Boudjema, B. Mele et al., Report of the Working Group on Standard Model 
Processes, in G. Altarelli, T. Sjostrand, F. Zwirner, eds.. Proceedings of the 
Workshop on Physics at LEP 2, CERN Yellow Report CERN 96-01 (1996), 
Vol. 1, 207-248, and references therein. 

[3] D. Bardin, M. Bilenky, D. Lehner, A. Olchevski and T. Riemann, Semi- 
analytical approach to four-fermion production in e~^e~ annihilation, in: 
T. Riemann and J. Bliimlein (eds.), Proc. of the Zeuthen Workshop on 
Elementary Particle Theory "Physics at LEP200 and Beyond", Nucl. Phys. 
B (Proc. Suppl.) 37B (1994), 148-157. 

[4] F.A. Berends, R. Pittau and R. Kleiss, Nucl. Phys. B424 (1994) 308. 

[5] D. Bardin and T. Riemann, Nucl. Phys. B462 (1996) 3. 

[6] D. Lehner, Ph.D. thesis, Humboldt-Universitat zu Berlin, Germany (1995), 
Internal Report DESY-IfH 95-07, e-Print Archive: |hep-ph/9512"301; 



D. Bardin, D. Lehner and T. Riemann, Nucl. Phys. B477 (1996) 27. 

[7] T. Muta, R. Najima and S. Wakaizumi, Mod. Phys. Lett. Al (1986) 203. 

[8] D. Bardin, A. Leike and T. Riemann, Phys. Lett. B353 (1995) 513. 

[9] D. Bardin, A. Leike and T. Riemann, Phys. Lett. B344 (1995) 383. 

[10] D. Bardin, M. Bilenky, A. Olchevski and T. Riemann, Phys. Lett. B 308, (1993) 
403; E: ibid, B 357, (1995) 725. 

[11] D. Bardin, W. Beenakker and A. Denner, Phys. Lett. B317 (1993) 213. 

[12] V.S. Fadin, V.A. Khoze, A.D. Martin and A. Stirling, Phys. Lett. B363 (1995) 
112. 

[13] Particle Data Group, R.M. Barnett et al., "Review of Particle Physics", Phys. 
Rev. D54 (1996). 

[14] V.S. Fadin, V.A. Khoze and A.D. Martin, Phys. Lett. B311 (1993) 311. 

[15] F.A. Berends, W. Beenakker et al., Report of the Working Group on WW cross- 
sections and distributions, in G. Altarelli, T. Sjostrand, F. Zwirner, eds., 
Proceedings of the Workshop on Physics at LEP 2, CERN Yellow Report CERN 
96-01 (1996), Vol. 1, 79-139, and references therein. 

[16] G. Gounaris, J.-L. Kneur and D. Zeppenfeldet al.. Report of the Working Group 
on Triple Gauge Boson Couplings, in G. Altarelli, T. Sjostrand, F. Zwirner, 
eds.. Proceedings of the Workshop on Physics at LEP 2, CERN Yellow Report 
CERN 96-01 (1996), Vol. 1, 525-576, and references therein. 



37 



[17] J. Biebel and T. Riemann, in preparation. 



[18] D. Bardin, M. Bilenky, A. Olchevski and T. Riemann, Fortran code 
GENTLE_1993, unpublished; 

D. Bardin, M. Bilenky, D. Lehner, A. Leike, A. Olchevski and T. Riemann, 
Fortran code GENTLE/4f an_1995; see also [1], p. 26, and [24], p.304, and D. 
Bardin, D. Lehner, A. Leike and T. Riemann, Internal Report DESY-Zeuthen 
96-05 (1996) [lhep-ph/9603438i . 

[19] A. Leike, preprint LMU-06-96, to appear in: T. Riemann and J. Bliimlein, 
eds., Proc. of the Zeuthen Workshop on Elementary Particle Theory "QED 
and QCD in higher orders", Rheinsberg, Germany, April 1996, e-Print Archive: 
|hep-ph/9606262| . 

[20] K. Hagiwara, S. Ishihara, R. Szalapski and D. Zeppenfeld, Phys. Lett. B283 
(1992) 353 and Phys. Rev. D48 (1993) 2182. 

[21] S. Eidelman and F. Jegerlehner, Z. Phys. C67 (1995) 585. 

[22] K. Hagiwara, K. Hikasa, R. Peccei and D. Zeppenfeld, Nucl. Phys. B282 (1987) 
253. 

[23] A. Pankov and N. Paver, ICTP Tileste preprint IC/96/175 [ |hep-ph/9610509i . 

[24] M. Mangano, R. Ridolfi et al.. Report of the Working Group on Event 
Generators for Discovery Physics, in G. Altarelli, T. Sjostrand, F. Zwirner, 
eds.. Proceedings of the Workshop on Physics at LEP 2, CERN Yellow Report 
CERN 96-01 (1996), Vol. 2, 299-353. 



38 



